arXiv: 1504.00271v2 [cond-mat.mtrl-sci] 2 Apr 2015 


Optical properties of inhomogeneous metallic hydrogen plasmas 

N. Van den Broeck,^ F. Brosens,^ J. Tempered and I.F. Silvera^ 

^ Theory of Quantum and Complex systems, Universiteit Antwerpen, Belgium 
^ Lyman Laboratory of Physics, Harvard University, Cambridge, Massachusetts 02138 

(Dated: April 3, 2015) 

We investigate the optical properties of hydrogen as it undergoes a transition from the insulating 
molecular to the metallic atomic phase, when heated by a pulsed laser at megabar pressures in 
a diamond anvil cell. Most current experiments attempt to observe this transition by detecting a 
change in the optical reflectance and/or transmittance. Theoretical models for this change are based 
on the dielectric function calculated for bulk, homogeneous slabs of material. Experimentally, one 
expects a hydrogen plasma density that varies on a length scale not substantially smaller than the 
wave length of the probing light. We show that taking this inhomogeneity into account can lead 
to significant corrections in the reflectance and transmittance. We present a technique to calculate 
the optical properties of systems with a smoothly varying density of charge carriers, determine the 
optical response for metallic hydrogen in the diamond anvil cell experiment and contrast this with 
the standard results. Analyzing recent experimental results we obtain = (2.1 ± 1.3) x 10® (fl 

cm)“® for the conductivity of metallic hydrogen at 170 GPa and 1250 K. 


INTRODUCTION 


Metallic hydrogen (MH) has been a challenge for ex¬ 
periment and theory since the initial proposal by Wigner 
and Huntington that under sufficient pressure hydrogen 
would dissociate to an atomic metallic solid [ij. In the 
past decades there have been a number of reports of pro¬ 
ducing metallic hydrogen in the laboratory (l-[l^. but 
only a shock experiment [1, and a recent static pres¬ 
sure experiment [I^ have provided convincing evidence 
of metallization. Many of the experiments use optical 
properties (when possible) to determine the state of the 
hydrogen, as this is a non-intrusive technique. In experi¬ 
ments where transmittance and reflectance are measured, 
the interpretation of optical properties is based on treat¬ 
ing the sample as a uniform metallic slab. Also current 
theoretical predictions ll|, 121 of the reflectance of MH 
are derived for a block of uniform density. 


In reality the samples are inhomogeneous conductors 
because the hot part of the sample is bounded by low 
temperature non-metallic regions of the same material. 
Here we show that the slab model introduces inaccuracies 
in the analyses of optical properties for an inhomogeneous 
plasma density and we propose a theoretical description 
that remedies this. Wave propagation through inhomo¬ 
geneous plasmas has been widely studied, but here we 
no longer assume that the inhomogeneity is on a length 
scale longer than the wavelength, and focus on the oppo¬ 
site case. We introduce an algorithm that also allows the 
presence of other materials (e.g. a tungsten absorber or a 
cladding layer) to be taken into account and that allows 
the effect of edge smoothing or diffusion of materials into 
each other to be investigated. The model is applicable to 
any situation where inhomogeneous plasma densities can 
occur, and thus is also relevant in, for example, studies 
of fusion plasmas, plasma materials processing reactors, 
and the ionosphere. Combining our method with the re¬ 


cent measurements of Zaghoo, Salamat and Silvera in a 
pulsed-laser heated diamond anvil cell (Ref. [l^, here¬ 
after referred to as ZSS), we estimate the conductivity 
of MH, a quantity that is of importance in planetary 
science as high-pressure high-temperature hydrogen is a 
major component of gas giants. 


THEORETICAL DESCRIPTION 


We consider a sample with charge carriers (here elec¬ 
trons) with a density profile n{r). When monochromatic 
light of frequency oj illuminates the sample, the total elec¬ 
tric field E(r,a;) causes the electrons to be displaced out 
of equilibrium over a distance X(r, w) that must be cal¬ 
culated in general from a microscopic response theory. 
Here, since the subject of our paper is not the micro¬ 
scopic theory but rather the macroscopic inhomogeneity 
in n(r), we will use a generalisation of the Drude response 
given in SI units by 


X(r,a;) 


e E(r, w) 

Ule 0j{uj -I- z/r) ’ 


( 1 ) 


where e and mg are the electron charge and mass, re¬ 
spectively, and T is the Drude relaxation time. The dis¬ 
placements X(r,a;) lead to induced charge and current 
densities given by 


uj) = -eV • [n(r)X(r, w)], (2) 

J*"'^(r, w) = iwen(r)X(r, w). (3) 


It is useful to rescale the density profile n(r) by the bulk 
free electron density uq, introducing /(r) = n{r)/no. 
The induced charge density can then be rewritten using 
the definition of relative permittivity of a Drude metal, 
e(a;) = + where Wpj = \Jme^ j{meSvac) 

is the plasma frequency, no the bulk free electron density. 
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and Cyac the vacuum permittivity: 

p™'^(r,a;) = [1 - £^( 0 ;)] V • [/j(r)E(r, w)] . (4) 

3 

Here we also introduced an index j in order to take into 
account the presence of different materials. The profile 
function fj (r) defines the location of material j and can 
describe smoothing of the edges of the material. Simi¬ 
larly, the induced current density m can be rewritten 
as 

J“'^(r,a;) = ^ [1 - ej(w)] /j(r)E(r, w). (5) 
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These induced charges and currents in turn induce elec¬ 
tric and magnetic fields. Note that in the above expres¬ 
sions, E(r, oj) is the total field, i.e. the sum of (i) the field 
induced by J™'^, and (ii) the externally applied field 
(the light beam). Relations (HI) and (O can be used to 
eliminate the induced charge and current density from 
the Maxwell equations. Writing these in terms of the 
vector potential A(r,w) and the scalar potential ^!)(r,w) 
instead of electric and magnetic field, we obtain the fol¬ 
lowing set of equations: 

^ (1 - £,) [fjA^ - yfj ■ + tcoA )], (6) 

j 

and 

Ta + ^ (1 - e,) /, {Vcj3 - tu3A) 

(7) 

with c the velocity of light in vacuum. Here and 
j^md represent the scalar and vector potentials induced 
by and whereas (j) and A represent the total 
scalar and vector potentials, including also the externally 
applied potentials. Note that we have left out the r, w de¬ 
pendence for notational simplicity. The Coulomb gauge 
V • A = 0 has been used. 

In a diamond anvil cell (DAC), equations (©-(iTl) are 
simplified by taking the uniaxial geometry of the system 
into account. Here, we take the x-axis to be aligned along 
the optical axis of the diamond anvils, and assume ho¬ 
mogeneity in the yz plane. The incoming light is consid¬ 
ered to be a plane wave travelling along the x-axis, cor¬ 
responding to = 0 and A®“* {x,u3) = 

When using this in the differential equations (O and ([T]) , 
it is readily seen that = 0. The only 

non-zero component A™'^ can be found from the remain¬ 
ing differential equation: 




A 


ind 

y 
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This equation is solved numerically on a grid of N points 
with spatial divisions Sx (a non-uniform grid is also pos¬ 
sible). Using the central finite difference method, the 
differential equation can be rewritten as the recursion 
relation 


Ua+l — X /1 (1 C) /j i^a) 


^ikxo 


k^Sx^ 


2 - (^1 - (1 - Cj)/i (a:^)) Ua - Ua-i (9) 


where Ua indicates the solution for A™‘^/Ao on grid point 
a. When the values of the two initial grid points ui 
and U 2 are known, this recursion relation can be used 
to calculate the remaining points. If we take the incom¬ 
ing wave to start from Xi and propagate towards Xn, 
then the values of the induced solution at the first two 
points will be the reflected wave. This statement as¬ 
sumes that the first grid points are located in the sur¬ 
rounding medium outside the system under considera¬ 
tion. In this case ui = Cr and U 2 = where 

Cr is an undefined complex constant connected to the 
reflectance through R = \Cr\^. Similarly, we define Ct 
from the last grid point, again located in the surround¬ 
ing medium, un{Cr) = , and relate it to the 

transmittance coefficient through T = \1 -I-CtI^, where 
the term 1 comes from the external wave. Furthermore, 
we must demand that the solution at un and un-i is in 
fact an outgoing wave of the form Hence we must 

impose un-i {Cr) = un{Cr). The goal of the 

algorithm will thus be to minimize 

\un-i (Cr) — UN {Cr)\ (10) 


with respect to Cr. The solution yields the correct Cr 
and Ct from which the reflectance R and transmittance 
T can be calculated. In turn, these two variables easily 
provide the absorption A = \ — R—T and thus completely 
define the optical response of the system. 


EFFECT OF INHOMOGENEITY 


When considering a semi-infinite block of material with 
sharp interface to the vacuum and permittivity e(a;), 
the reflectance is given by the textbook result Ro{uj) = 
K-y/e — l)/( y ^ -I- l)p . This formula has been used in 
Refs. [11LI12|| in conjunction with microscopic derivations 
of the permittivity e(a;) of metallic hydrogen to predict 
the change in reflectance upon metallization. In order 
to demonstrate that Ro{uj) can give an incorrect result, 
and the formalism presented above provides a different 
result appropriate for an inhomogeneous density profile, 
we compute the reflectance for a density profile modeled 
by 


fix) 


1 _ exp {£) _^ 

2 cosh (^) -I- cosh [(2a; — f)/«:] 


( 8 ) 


( 11 ) 
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Figure 1. (Color online) The reflectance of the density profile 
presented in eq. GD) for MH in function of the thickness d for 
K = 0, 5, 10 and 15 nm (disks, squars, triangles and diamonds 
resp.). The solid line coinciding with the disks represents the 
analytical result for a hard wall. The black dashed line shows 
the analytical transmittance. Transmittance for different k is 
not shown since they all coincide with the dashed line. Data 
points with combinations of d and k resulting in a density 
prohle that does not reach 1 have been omitted. Inset (a) 
shows the reflectance as a function of k for d = 500 nm. In 
this case, transmittance is zero for all ft. Inset (b) shows 
the density profile for d = 60, 100 and 140 nm with ft = 5 
nm. The parameters used in all these hgures are: permittivity 
e = —1.47 + 13. 6i and A = 500 nm. 


This represents a layer of thickness d, with edges 
smoothed over a distance ft, as illustrated in inset (b) 
of Fig. [TJ In the limit ft —> 0 we retrieve a sharp-edged 
layer. In Fig. [T]the reflectance for such a smooth-edged 
layer is shown as a function of increasing thickness d, 
keeping ft fixed (at 0, 5, 10 and 15 nm for disks, squares, 
triangles and diamonds respectively). As the layer thick¬ 
ness grows, the reflectance goes to a constant value, but 
only for ft —f 0 does this correspond to the result 
for a sharp interface, shown as a full curve. As the in¬ 
terface is smoothed, the reflectance decreases in favor 
of the absorptance. This decrease is also shown in in¬ 
set (a). This means that even for thick films, the re¬ 
flectance is influenced strongly by the smoothness of the 
interface. The enhancement of absorption can be related 
to a smooth changing permittivity, which is comparable 
to techniques used in graded refractive index antireflec¬ 
tion coatings [l^ . 

In Fig. [21 results are shown for a different kind of den¬ 
sity profile. In pulsed-laser heated DACs one expects the 
temperature to drop off from a maximum at the heater 
to zero smoothly 1^, over a distance 5. A similar pro¬ 
file can then be expected for the density of dissociated 
molecules. This type of profile, illustrated in the inset of 
Fig. 121 and parametrized by f{x) = exp[l — exp(—cc/J)] 
for a: < 0, is used to plot the reflectance and transmit¬ 


Figure 2. (Color online) Reflectance and transmittance are 
shown as a function of layer thickness for a metal with 
e = —1.47 -I- 13.6i. The solid (reflectance) and dashed (trans¬ 
mittance) curves show the result for a homogeneous slab of 
thickness d, and the disks show the result for the inhomoge¬ 
neous density profile f(x) shown in the inset and representing 
a smooth increase from zero to maximum density over a dis¬ 
tance 5, after which the density is zero (i.e. for a: > 0). Light 
with wavelength A = 500 nm is incident from x = —oo. 


tance as a function of (5 in Fig. |21 The disks and squares 
show the results with the current formalism, and the solid 
and dashed curves show the analytical result for a sharp- 
edged layer with the same total amount of charge car¬ 
riers. For thin films, the precise shape of the density 
profile is unimportant since the layer becomes effectively 
a two-dimensional sheet of charge carriers. However, as 
the thickness is increased, we again find that the inho¬ 
mogeneity of the density profile leads to a reduction of 
reflectance in favor of absorption. Hence, it is necessary 
to take the inhomogeneity into account when inferring 
the him thickness from an optical rehectance measure¬ 
ment, once the him becomes thicker than a few tens of 
nanometer. 


METALLIC HYDROGEN IN A DAC 


In ZSS the phase transition to liquid metallic hydrogen 
was observed in a DAC by heating the sample along an 
isochore. A thin him of tungsten is used to heat the 
hydrogen by absorbing a few hundreds nanosecond long 
laser pulse 1^. The W him has been constructed to 
absorb 50% of the incoming light and transmit another 
25%, which corresponds to a him thickness of 8 nm for a 
laser of A = 500 nm. We take the frequency dependent 
tungsten permittivity from tables published in [l5| . The 
tungsten him is separated from the hydrogen by a 5 nm 
cladding layer of (optically inactive) alumina. 

As the pulsed laser power is increased, the tempera¬ 
ture of the tungsten him grows, and when the tempera- 
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ture T reaches the dissociation temperature T^, a layer 
of metallic hydrogen is formed, reducing transmittance 
and increasing reflectance from the sample. The den¬ 
sity of the metallic hydrogen is expected to be highest 
where it is in contact with the tungsten film. Combining 
the Arrhenius equation [l^ for the dissociation rate of 
molecules, and the temperature profile for pulsed laser 
heating derived in 0 , we propose the following density 
profile for the metallic hydrogen: 


fn (a;) = exp 




( 12 ) 


The maximum number density of metallic hydrogen is 
taken to be twice the molecular hydrogen density ob¬ 
tained from the equation of state [l7| . The maximum 
MH density occurs at a; = 0, close to the tungsten layer, 
and drops off for a; > 0. In the region a; < 0, the density 
of metallic hydrogen is zero, and for the remainder of the 
calculations we set T = T^. The resulting density profile 
(fT^ is similar to those illustrated in the inset of Fig. [2] 

The theoretical description introduced above enables 
us to investigate the difference in reflectance (transmit¬ 
tance) between this inhomogeneous density profile and 
the reflectance (transmittance) obtained for a uniform 
slab of MH. This is of importance in order to correctly 
derive the properties of the newly detected phase from 
the observed changes in optical response. From Ref. 111 
we have r = 4.8 x 10“^^ s and ujpi = 22 eV at 30000 
K and 152 GPa, resulting in no = 3.5 x 10^® m“^ and 
a permittivity cmh = —1.47 -I- 13.6f at A = 500 nm as 
modelled by e{uj) = 1 — Wpj/[w(a; -I- i/r)]. 

Figure [3] contrasts the results for a uniform layer of hy¬ 
drogen with that of a density profile as given by eq. (113. 
The change in reflectance (transmittance) is calculated as 
the ratio of the reflectance (transmittance) with a hydro¬ 
gen layer of thickness 5th and without the hydrogen layer 
{5th 0). Results for a uniform hydrogen density are 
shown as a solid curve (dashed curve), and results for 
the non-uniform case are shown as disks (squares). The 
non-uniform layer of the thickness 5th is compared to a 
uniform him with thickness d chosen such that the den¬ 
sity integrated over x is the same. Discrepancies start 
appearing already at ca. 10 nm thickness, and are more 
pronounced for the change in rehection. This is due to 
the fact that the inhomogeneous density prohle leads to a 
stronger absorption, which may even result in a relative 
decrease of the rehectance for thicker layers. In prac¬ 
tice, this increased absorption by the metallic hydrogen 
layer will lead to an increase in temperature above the 
dissociation temperature. 

We can use these results to analyse the optical mea¬ 
surement of ZSS. We use the data at 170 GPa, namely 
the reflectance change at 980 nm and the transmittance 
change at 980 nm and 633 nm, and focus on data in 
the temperature range 1260-1290 K just above the heat¬ 
ing plateau that indicates a first order transition to the 
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Figure 3. (Color online) The change in reflectance and trans¬ 
mittance at A = 500 nm, is shown as a function of the thick¬ 
ness d of a layer of metallic hydrogen on top of a 15 nm thick 
layer of tungsten. The disks and squares show the reflectance 
and transmittance, respectively, for an inhomogeneous den¬ 
sity profile as given by eq. dnj. The solid and dashed line 
show reflectance and transmittance, respectively, for a layer 
of hydrogen with uniform density. The thickness of the uni¬ 
form layer is chosen such that the total amount of electrons is 
the same as that for the non-uniform profile, and it is shown 
on the upper axis. 


plasma phase 0. We set ujpi = 26 eV, corresponding 
to a charge density twice the hydrogen molecular density 
at 170 GPa, and use a fit to the observed reflectance 
and transmittances to obtain a Drude relaxation rate 
T = (1.1±0.4) X 10“^^ sec. The corresponding Drude con¬ 
ductivity is = (2100±1300) (D cm)“^. This value 

is in agreement with the value of 2000 (D cml“^ found for 
fluid hydrogen in shock wave experiments [9| at 140 GPa 
and 2600 K, even though the shock experiment finds a 
continuous transition from a semiconducting to a metallic 
fluid, whereas ZSS observe a first-order phase transition. 
This may indicate that in both cases the metallic fluid 
is very disordered, and has a conductivity close to the 
Mott-Ioffe-Regel minimum metal conductivity[l^. The 
relaxation rate corresponding to the minimum metallic 
conductivity in two dimensions is T^in = 0.68 x 10“^^ 
sec. The film thickness compatible with our fit is 1-3 
nm, which also hints at the 2D nature of the MH film 
formed in the experiment. Finally, note that changing 
the plasma frequency to 21 eV (the value reported in 
01 ) would yield t = (1.7 ± 0.9) x 10“^^ sec. 

CONCLUSIONS 

In this article we propose a method to obtain the re¬ 
flectance, transmittance and absorptance for an inhomo¬ 
geneous charge density, given the microscopic response 
(i.e. the dielectric function). The results are contrasted 
to results for uniform slabs of materials obtained with 
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the standard Maxwell boundary conditions for sharp in¬ 
terfaces. The theoretical description is then applied to 
the case of liquid metallic hydrogen [l^ in a DAC, where 
it is found that taking this into account reduces the re¬ 
flectance of the metallic hydrogen with respect to that 


obtained in previous works for thick samples [11 


and increases the absorptance. Finally, we extract the 
conductivity of metallic hydrogen from the measurement 
data of ZSS. The theory presented here will be important 
for forthcoming experiments that attempt to refine the 
optical measurements used to probe the newly detected 
metallic plasma phase of hydrogen. 
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